arXiv:1509.04670v2 [hep-ex] 11 Nov 2015 


EUROPEAN ORGANISATION FOR NUCLEAR RESEARCH (CERN) 




Phys. Rev. D92, 092004 (2015) 
DOI: 10.1103/PhysRevD.92.092004 


CERN-PH-EP-2015-225 
12th November 2015 


Searches for Higgs boson pair production in the 
hh bbTT, yyWW*, yybb, bbbb channels with the ATLAS 


detector 


The ATLAS Collaboration 


Abstract 


Searches for both resonant and nonresonant Higgs boson pair production are performed in 
the hh —> bbTT, yyWW* final states using 20.3 fb”^ of pp collision data at a center-of-mass 
energy of 8 TeV recorded with the ATLAS detector at the Large Hadron Collider. No evid¬ 
ence of their production is observed and 95% confidence-level upper limits on the produc¬ 
tion cross sections are set. These results are then combined with the published results of the 
hh^yybb, bbbb analyses. An upper limit of 0.69 (0.47) pb on the nonresonant hh produc¬ 
tion is observed (expected), corresponding to 70 (48) times the SM gg —> hh cross section. 
For production via narrow resonances, cross-section limits of hh production from a heavy 
Higgs boson decay are set as a function of tbe heavy Higgs boson mass. The observed (ex¬ 
pected) limits range from 2.1 (1.1) pb at 260 GeV to 0.011 (0.018) pb at 1000 GeV. These 
results are interpreted in the context of two simplified scenarios of the Minimal Supersym¬ 
metric Standard Model. 
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1 Introduction 


The Higgs boson discovered at the LHC in 2012 [1, 2] opens a window for testing the scalar sector of 
the Standard Model (SM) and its possible extensions. Since the discovery, significant progress has been 
made in measuring its coupling strengths to fermions and vector bosons [3-6] as well as in studying 
its spin and its charge-conjugate and parity (CP) properties [7, 8]. All results are consistent with those 
expected for the SM Higgs boson (here denoted by h). Within the SM, the existence of the Higgs boson is 
a consequence of the electroweak symmetry breaking (EWSB). This also predicts self-coupling between 
Higgs bosons, the measurement of which is crucial in testing the mechanism of EWSB. The self-coupling 
is one mechanism for Higgs boson pair production as shown in Eig. 1(a). Higgs boson pairs can also be 
produced through other interactions such as the Higgs-fermion Yukawa interactions (Eig. 1(b)) in the 
Standard Model. These processes are collectively referred to as nonresonant production in this paper. 



Figure 1; Leading-order Feynman diagrams of the nonresonant production of Higgs boson pairs in the Standard 
Model through (a) the Higgs boson self-coupling and (b) the Higgs-fermion Yukawa interactions. 

Higgs boson pair production at the LHC as a probe of the self-coupling has been extensively studied in 
the literature [9-13]. One conclusion [14] is that the data collected so far (approximately 25 fb“^ in total) 
are insensitive to the self-coupling in the SM, because of the expected small signal rates [15-17] and large 
backgrounds. However, it is essential to quantify the sensitivity of the current dataset and to develop tools 
for future measurements. Moreover, physics beyond the Standard Model (BSM) can potentially enhance 
the production rate and alter the event kinematics. For example, in the Minimal Supersymmetric Standard 
Model (MSSM) [ 1 8] , a heavy CP-even neutral Higgs boson H can decay to a pair of lighter Higgs bosons. 
Production of H followed by its decay H^hh would lead to a new resonant process of Higgs boson pair 
production, in contrast to the nonresonant hh production predicted by the SM (Fig. 1). In composite Higgs 
models such as those discussed in Refs. [19, 20], increased production of nonresonant Higgs boson pairs 
is also expected. 

Both the ATEAS and CMS collaborations have searched for nonresonant and/or resonant Higgs boson 
pair production [21-23]. In particular, ATEAS has published the results of searches in the hh^yybb [21] 
and hh —> bbbb [22] decay channels.' In this paper, searches in two additional hh decay final states, 
bbTT and yyWW*, are reported. For the hh bbTT analysis, one tau lepton is required to decay to an 
electron or a muon, collectively referred to as and the other tau lepton decays to hadrons (Thad)- For 
hh yyWW*, the h —> WW* —> £vqq' decay signature is considered in this study. The results of these 
new analyses are combined with the published results of hh —> yybb and hh bbbb for both nonresonant 
and resonant production. The resonance mass mu considered in this paper ranges from 260 GeV to 


* Notations indicating particle charges or antiparticles are generally omitted throughout this paper. 
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1000 GeV. The lower bound is dictated by the 2m/, threshold while the upper bound is set by the search 
range of the hh bbxT analysis. The light Higgs boson mass m/, is assumed to be 125.4 GeV, the 
central value of the ATLAS measurement [24]. At this mass value, the SM predictions [25-27] for the 
decay fractions of hh^bbbb, bbTT, bbyy and yyWW* are, respectively, 32.6%, 7.1%, 0.26% and 0.10%. 
The resonant search assumes that gluon fusion is the production mechanism for a heavy Higgs boson 
that can subsequently decay to a pair of lighter Higgs bosons, i.e., gg ^ H ^ hh. Furthermore, the 
heavy Higgs boson is assumed to have a width significantly smaller than the detector resolution, which 
is approximately 1.5% in the best case (the hh —> yybb analysis). The potential interference between 
nonresonant and resonant production is ignored. 

This paper is organized as follows. For the hh bbTT and hh yy'W'W* analyses, data and Monte 
Carlo (MC) samples are described in Sec. 2 and the object reconstruction and identification are outlined 
in Sec. 3. In Secs. 4 and 5, the separately published hh —> yybb and hh —> hhhh analyses are briefly sum¬ 
marized. The hh bbTT and hh yy'W'W* analyses including even! selecfion, background estimations 
and systematic uncertainties are presented in Secs. 6 and 7, respectively. The statistical and combination 
procedure is described in Sec. 8. The results of the hh bbTT and hh —> yyWW* analyses, as well as their 
combinations with the published analyses are reported in Sec. 9. The implications of the resonant search 
for two specific scenarios of fhe MSSM, hMSSM [28, 29] and low-fb-high [30], are discussed in Sec. 10. 
These scenarios make specific assumpfions and/or choices of MSSM paramefers fo accommodafe fhe 
observed Higgs boson. Finally, a summary is given in Sec. 11. 


2 Data and Monte Carlo samples 

The dafa used in fhe searches were recorded in 2012 wifh fhe ATLAS defecfor af fhe Large Hadron 
Collider in profon-profon collisions af a cenler-of-mass energy of 8 TeV and correspond fo an infegrafed 
luminosity of 20.3 fb“^ The ATLAS defector is described in defail in Ref. [31]. Only dafa recorded when 
all subdefecfor systems were properly functional are used. 

Signal and background MC samples are simulafed wifh various even! generafors, each inferfaced fo Py- 
THiA v8.175 [32] for parfon showers, hadronizafion and underlying-evenf simulation. Parfon disfribufion 
functions (PDFs) CTIO [33] or CTEQ6L1 [34] for fhe profon are used depending on fhe generator in 
quesfion. MSTW2008 [35] and NNPDF [36] PDFs are used fo evaluafe sysfemafic uncerfainfies. Table 1 
gives a brief overview of fhe evenf generafors, PDFs and cross secfions used for fhe hh —> bbTT and 
hh —> yyVTVT* analyses. All MC samples are passed fhrough fhe ATLAS defecfor simulafion program [37] 
based on GEANT4 [38]. 

Signal samples for bofh nonresonanf and resonanf Higgs boson pair producfion are generafed using fhe 
leading-order MadGraph5 vL5.14 [39] program. The nonresonanf producfion is modeled using fhe SM 
DiHiggs model [40, 41] while fhe resonanf producfion is realized using fhe HeavyScalar model [42], 
bofh implemented in MadGraph5. The heavy scalar H is assumed fo have a narrow decay widfh of 
10 MeV, much smaller fhan fhe experimenfal resolution. The SM predicfion for fhe nonresonanf gg^hh 
producfion cross section is 9.9 ± 1.3 fb [17] wifh m/, = 125.4 GeV from fhe nexl-fo-nexl-lo-leading-order 
calculation in QCD. 

Single SM Higgs boson producfion is considered as a background. The Powheg rl655 generator [43- 
45] is used fo produce gluon fusion (ggP) and vector-boson fusion (VBE) evenfs. This generator calcu¬ 
lates QCD correcfions up fo nexf-fo-leading order (NEC), including finife bottom- and fop-quark mass 
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Table 1: List of MC generators and parton distribution functions of the signal and background processes used by 
the hh —> bbTT and hh yyWW* analyses. SM cross sections used for the normalization are also given. For the VTZ 
and ZZ processes, contributions from y* are included and the cross sections quoted are for mziy> > 20 GeV. 


Process 

Event generator 

PDL set 

Cross section [pb] 

Background processes 

P-i-jets 

Alpgen -i- PythiaS 

CTEQ6L1 

normalized to data 

Diboson: WW 

Powheg -i- PythiaS 

CTIO 

55.4 

Diboson: WZ 

Powheg -i- PythiaS 

CTIO 

22.3 

Diboson: ZZ 

Powheg -i- PythiaS 

CTIO 

7.3 

tt 

Powheg -i- PythiaS 

CTIO 

253 

Single top: t-channel 

AcerMC -i- PythiaS 

CTEQ6L1 

87.8 

Single top: ^-channel 

Powheg -i- PythiaS 

CTIO 

5.6 

Single top: Wt 

Powheg -i- PythiaS 

CTIO 

22.0 

gg^h 

Powheg -i- PythiaS 

CTIO 

19.2 

qq' ^qq'h 

Powheg -i- PythiaS 

CTIO 

1.6 

qq^Vh 

PythiaS 

CTEQ6L1 

1.1 

qqlgg^tth 

PythiaS 

CTEQ6L1 

0.13 

Signal processes 

Nonresonant gg —> hh 

MadGraph5 -I- PythiaS 

CTEQ6L1 

0.0099 

Resonant gg^H^hh 

MadGraph5 -I- PythiaS 

CTEQ6L1 

model dependent 


effects [46]. The Higgs boson transverse momentum (px) spectrum of the ggF process is matched to 
the calculated spectrum at next-to-next-to-leading order (NNLO) and next-to-next-to-leading logarithm 
(NNLL) [47] in QCD corrections. Events of associated production qq Vh (here V = W or Z) and 
qqigg —> tth are produced using the PythiaS generator [32]. All of these backgrounds are normalized using 
the state-of-the-art theoretical cross sections (see Table 1) and their uncertainties compiled in Refs. [25- 
27]. 

The Alpgen v2.1.4 program [48] is used to produce the P-f-jets samples. The Powheg generator is used to 
simulate top quark pair production {tt) as well as the ^-channel and Wt processes of single top production; 
the f-channel process of single top production is simulated using the AcerMC v38 program [49]. The tt 
cross section has been calculated up to NNLO and NNLL in QCD corrections [50]. Cross sections for the 
three single-top processes have been calculated at (approximate) NNLO accuracy [51-53]. The Powheg 
generator is used to simulate diboson backgrounds {WW, fPZ and ZZ). The diboson production cross 
sections are calculated at NLO in QCD corrections using the MCLM program [54, 55]. 


3 Object identification 


In this section, object reconstruction and identification for the hh —> bbTT and hh —> yyWW* analyses 
are discussed. The hh bbTT and hh yy'W'W* analyses are developed following the h ^ tt [6] and 
h^yy [56] studies of single Higgs bosons, respectively, and use much of their methodology. 


4 











Electrons are reconstructed from energy clusters in the electromagnetic calorimeter matched to tracks in 
the inner tracker. The calorimeter shower profiles of electron candidates must be consistent with those 
expected from electromagnetic interactions. Electron candidates are identified using fighf and medium 
criferia [57] for fhe hh —> bbTT and hh yyWW* analyses, respectively. The selecfed candidafes are 
required fo have fransverse momenfa^ pj > 15 GeV and be wifhin fhe defecfor fiducial volume of \r]\ < 
2.47 excluding 1.37 < \r]\ < 1.52, fhe fransifion region befween fhe barrel and endcap calorimeters. 
Muons are idenlified by mafching tracks or segments reconstructed in the muon spectrometer with tracks 
reconstructed in the inner tracker. They are required to have px > 10 GeV and \t]\ < 2.5. Both the 
electrons and muons must satisfy calorimeter and track isolation requirements. The calorimeter isolation 
requires that the energy deposited in the calorimeter in a cone of size AR = = 0.2 around 

the lepton (electron or muon), excluding the energy deposited by the lepton itself, is less than 6% (20%) 
of the px of the lepton for the hh^bbzT (hh^yyWW*) analysis. The track isolation is defined similarly: 
fhe scalar px sum of additional fracks originating from fhe primary verfex wifh px > 1 GeV in a cone of 
size AR = 0.4 around fhe lepfon is required fo be less fhan 6% (15%) of fhe px of fhe lepton frack for fhe 
hh bbTT {hh yyWW*) analysis. 

Photons are reconsfrucfed from energy clusfers in fhe elecfromagnefic calorimefer wifh fheir shower pro¬ 
files consisfenf wifh elecfromagnefic showers. A significanf fracfion of phofons converf info e'*'e~ pairs 
inside fhe inner fracker. Thus photon candidafes are divided info unconverted and converted categories. 
Clusfers wifhouf mafching fracks are considered as unconverted phofons, while clusters mafched fo fracks 
consisfenf wifh conversions are considered as converfed phofons. Phofon candidafes musf fulfill fhe fighf 
idenfificafion criferia [58] and are required fo have px > 25 GeV and \r]\ < 2.37 (excluding fhe fransifion 
region 1.37 < Ipl < 1.52) and musf satisfy bofh calorimefer and frack isolation. The calorimeter isolation 
requires fhe additional energy in a cone of AR = 0.4 around fhe phofon candidafe fo be less fhan 6 GeV 
while fhe frack isolafion requires fhe scalar px sum of addifional fracks in a cone of AR = 0.2 around fhe 
phofon fo be less fhan 2.6 GeV. 

Jefs are reconsfrucfed using fhe anfi-k, algorifhm [59] wifh a radius parameter of = 0.4. Their energies 
are corrected for fhe average confribufions from pileup inferacfions. Jefs are required fo have px > 30 GeV 
and Ipl < 4.5. Eor fhe hh^yyWW* analysis, a lower px requiremenf of 25 GeV is applied for jefs in fhe 
cenfral region of Ipl < 2.4. To suppress confribufions from pileup inferacfions, jefs wifh px < 50 GeV and 
wifhin fhe accepfance of fhe inner fracker (Ipl < 2.4) musf have over 50% of fhe scalar sum of fhe px of 
fheir associated fracks confribufed by fhose originafing from fhe primary verfex. Jefs confaining f^-hadrons 
are idenfified using a mulfivariafe algorifhm (ft-fagging) [60]. The algorifhm combines information such 
as fhe explicif reconsfrucfion of fhe secondary decay vertices and frack impacf-paramefer significances. 
The operafing poinf chosen for bofh hh —> bbTT and hh —> yyWW* analyses has an efficiency of 80% for 
fhe Z>-quark jefs in ti evenfs. 

Hadronically decaying r candidafes (rhad) are reconsfrucfed using clusfers in fhe elecfromagnefic and 
hadronic calorimefers [61]. The fau candidafes are required fo have px > 20 GeV and \r]\ < 2.5. The 
number of fracks wifh px > 1 GeV associated wifh fhe candidafes musf be one or fhree and fhe fofal 
charge defermined from fhese fracks musf be +1. The fau idenfificafion uses calorimefer clusfer as well 
as fracking-based variables, combined using a boosfed-decision-lree mefhod [61]. Three working poinfs. 


^ ATLAS uses a right-hand coordinate system with the interaction point as its origin and the beam line as its z-axis. The x-axis 
points to the center of the LHC ring and y-axis points upwards. The pseudorapidity 77 is defined as 77 = - In tan(0/2), where 
8 is the polar angle measured with respect to the z-axis. The transverse momentum pj is calculated from the momentum p: 
Pt = p sin 8. 
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labeled loose, medium and tight [61], corresponding to different identification efficiencies are used. Ded¬ 
icated algorithms that suppress electrons and muons misidentified as Thad candidates are applied as well. 

The missing transverse momentum (with magnitude is the negative of the vector sum of the trans¬ 
verse momenta of all photon, electron, muon, Thad^ and jet candidates, as well as the px of all calori¬ 
meter clusters not associated with these reconstructed objects, called the soft-term contribution [62]. The 
hh^bbTT analysis uses the version of the calculation in the h^TT analysis [6]. In this calculation, 
the soft-term contribution is scaled by a vertex fraction, defined as fhe ratio of fhe summed scalar px of 
all fracks from fhe primary verfex nof mafched wifh fhe reconsfrucfed objecfs fo the summed scalar px 
of all tracks in the event. The hh yyWW* analysis, on the other hand, uses the £'™®®-significance em¬ 
ployed by the h^yy study [56]. It is defined as fhe ratio of fhe measured to its expected resolution 
estimated using the square root of the scalar sum of the transverse energies of all objects contributing to 
the calculation. 

4 Summary of hh yybb 

The hh —> yybb analysis, published in Ref. [21], largely follows the ATLAS analysis of the Higgs boson 
discovery in the h ^ yy decay channel [1, 56]. The search is performed in the ^/s = S TeV dataset 
corresponding to an integrated luminosity of 20.3 fb“^ The data were recorded with diphoton triggers 
that are nearly 100% efficient for events satisfying the photon requirements. Events must contain two 
isolated photons. The px for the leading (subleading) photon must be larger than 35% (25%) of the 
diphoton invariant mass niyy, which itself must be in the range of 105 < myy < 160 GeV. Events must 
also contain two energetic (^-tagged jets; the leading (subleading) jet must have px > 55 (35) GeV, and 
the dijet mass must fall within a window 95 < rntb < 135 GeV, as expected from the h ^ bb decay. A 
multivariate fj-tagging algorithm [60] that is 70% efficient for the fj-quark jets in ft events is applied. 

Backgrounds for both the resonant and nonresonant analyses are divided into two categories: events con¬ 
taining a single real Higgs boson (with h yy), and the continuum background of events not containing a 
Higgs boson. The former are evaluated purely from simulation, and are small compared to the continuum 
background, which is evaluated from data in the diphoton mass sidebands (the myy range of 105-160 GeV 
excluding the region of m/, ± 5 GeV). In the nonresonant analysis, an unbinned signal-plus-background fit 
is performed on the observed myy distribution, with the background from single Higgs bosons constrained 
to the expectation from the SM. The continuum background is modeled with an exponential function; the 
shape of the exponential function is taken from data containing a diphoton and dijet pair where fewer than 
two jets are f?-tagged. 

The resonant search proceeds in a similar manner, although it is ultimately a counting experiment, with an 
additional requirement on the four-object invariant mass myybb, calculated with the bb mass constrained to 
nih- This requirement on myybb varies with the resonance mass hypothesis under evaluation, and is defined 
as fhe smallesf window containing 95% of the signal events based on MC simulation. As in the nonres¬ 
onant search, the number of background events with real Higgs bosons is estimated from simulation. The 
continuum background in the myy signal window is extrapolated from the diphoton mass sidebands. A 
resonance with mass between 260 GeV and 500 GeV is considered in the search. 

The small number of events (nine) in the diphoton mass sideband leads to large statistical uncertainties 
(33%) on the dominant continuum background, so that most systematic uncertainties have a small effect 
on the final result. Eor the resonant search, however, systematic uncertainties with comparable effect 
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remain. Uncertainties of 0-30% (depending on the resonance mass hypothesis under consideration) are 
assigned due to the modeling of the myyhb shape using the data with less than two ( 7 -tagged jets. Additional 
uncertainties of 16-30% arise from the choice of functional form used to parameterize the shape of 

fltyybb- 

In the nonresonant analysis, extrapolating the sidebands into the diphoton mass window of the signal 
selection leads to a prediction of 1.3 continuum background events. An additional contribution of 0.2 
events is predicted from single Higgs boson production. A total of five events are observed, representing 
an excess of 2.4 standard deviations (cr). A 95% confidence level (CL) upper limit of 2.2 (1.0) pb is 
observed (expected) for cr{gg —> hh), the cross section of nonresonant Higgs boson pair production. For 
the resonant searches, the observed (expected) upper limits on crigg^H) x BR(H^hh) are 2.3 (1.7) pb 
at niH - 260 GeV and 0.7 (0.7) pb at mn - 500 GeV. 


5 Summary of hh bbbb 


The hh —> bbbb analysis [22] benefits from the large branching ratio of h^ bb. The analysis employs 
resolved as well as boosted Higgs boson reconstruction methods. The resolved method attempts to re¬ 
construct and identify separate f7-quark jets from the h ^ bb decay, while the boosted method uses a jet 
substructure technique to identify the h ^ bb decay reconstructed as a single jet. The latter is expected 
if the Higgs boson h has a high momentum. The boosted method is particularly suited to the search for 
a resonance with mass above approximately 1000 GeV decaying to a pair of SM Higgs bosons. For the 
combinations presented in this paper, resonances below this mass are considered and the resolved method 
is used as it is more sensitive. 

The analysis with the resolved method searches for two back-to-back and high-momentum bb systems 
with their masses consistent with m/, in a dataset at = 8 TeV corresponding to an integrated luminosity 
of 19.5 fb“' for the triggers used. The data were recorded with a combination of multijet triggers using 
information including the f 7 -quark jet tagging. The trigger is >99.5% efficient for signal events satisfying 
the offline selection. Candidate evenfs are required fo have af leasf four fj-lagged jefs, each wifh px > 
40 GeV. As in fhe hh —> yybb analysis, a mulfivariafe f?-lagging algorifhm [60] wifh an esfimafed efficiency 
of 70% is used fo fag jefs confaining f 7 -hadrons. The four highesf-px (^-fagged jefs are fhen used fo form 
fwo dijef sysfems, requiring fhe angular separation A/? in (p, (p) space befween fhe fwo jefs in each of fhe 
fwo dijef sysfems fo be smaller fhan 1.5. The fransverse momenfa of fhe leading and subleading dijef 
sysfems musf be greafer fhan 200 GeV and 150 GeV, respecfively. These selection criferia, driven parfly 
by fhe corresponding jef frigger fhresholds and parfly by fhe necessify fo suppress fhe backgrounds, lead 
fo significanf loss of signal accepfance for lower resonance masses. The resonanf search only considers 
masses above 500 GeV. The leading {mi 2 ) and subleading (^ 34 ) dijef invarianf mass values are required 
fo be consisfenf wifh fhose expecfed for fhe hh —> bbbb decay, salisfying fhe requiremenf: 


A 


r 0 

mi2 - mV, ) 


-I- 


^ 0 

m34 - m34 

, 0-34 


< 1 . 6 . 


Here (124 GeV) and m^^ (115 GeV) are fhe expecfed peak values from simulation for fhe leading 
and subleading dijef pair, respecfively, and cri 2 and 0^34 are fhe dijef mass resolufions, esfimafed from 
fhe simulation fo be 10% of fhe dijef mass values. More defails abouf fhe selecfion can be found in 
Ref. [22]. 
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After the full selection, more than 90% of the total background in the signal region is estimated to be 
multijet events, while the rest is mostly tt events. The Z+jets background constitutes less than 1% of the 
total background and is modeled using MC simulation. The multijet background is modeled using a fully 
data-driven approach in an independent control sample passing the same selection as the signal except 
that only one of the two selected dijets is Z>-tagged. This control sample is corrected for the kinematic 
differences arising from the additional f7-tagging requirements in the signal sample. The tt contribution is 
taken from MC simulations normalized to data in dedicated control samples. 

The dominant sources of systematic uncertainty in the analysis are the f^-tagging calibration and the 
modeling of the multijet background. The degradation in the analysis sensitivity from these uncertainties 
is below 10%. Other sources of systematic uncertainty include the tt modeling, and the jet energy scale 
and resolution, which are all at the percent level. A total of 87 events are observed in the data, in good 
agreement with the SM expectation of 87.0 ± 5.6 events. Good agreement is also observed in the four-jet 
invariant mass distribution, thus there is no evidence of Higgs boson pair production. For the nonresonant 
search, both the observed and expected 95% CL upper limit on the cross section cr{pp -^hh^ bbbb) is 
202 fb. For the resonant search, the invariant mass of the four jets is used as the final discriminanf from 
which fhe upper limif on fhe pofenfial signal cross secfion is exfracfed. The resulting observed (expected) 
95% CL upper limit on cr{pp^H^hh^bbbb) ranges from 52 (56) fb, at tuh = 500 GeV, to 3.6 (5.8) fb, 
at niH - 1000 GeV. 


6 hh-^bbTT 

This section describes the search for Higgs boson pair production in the hh bbrr decay channel, where 
only the final slate where one fau lepfon decays hadronically and fhe ofher decays lepfonically, bbbThud, 
is used. The dala were recorded wilh friggers requiring al leasl one lepfon wilh pj > 24 GeV. These 
Iriggers are nearly 100% efiicienl for evenfs passing fhe final selection. Candidate bbbThad evenfs are 
selected by requiring exacfly one lepfon wilh pj > 26 GeV, one hadronically decaying fau lepfon of fhe 
opposile charge wilh pj > 20 GeV meeling fhe medium criteria [61], and Iwo or more jels wilh px > 
30 GeV. In addifion, belween one and Ihree of fhe selected jels musf be fj-lagged using fhe multivariate 
b-taggei. The upper bound on fhe number of b-t&gged jels is designed fo make Ihis analysis sfalisfically 
independenl of fhe hh bbbb analysis summarized in Sec. 5. These criteria are colleclively referred fo 
as fhe “preselecfion”. 

The backgrounds from IT-i-jefs, Z —> tt, diboson {WW, fkZ and ZZ) and lop quark (bofh tt and single 
lop quark) produclion dominale fhe surviving sample and Iheir conlribufions are derived from a mixture 
of data-driven methods and simulation. The contribution from events with a jet misidentified as a Thad, 
referred to as the fake Thad background, are estimated using data with a “fake-factor” method. The method 
estimates contributions from VF-i-jets, multijet, Z-i-jets and top quark events that pass the event selection 
due to misidentified Thad candidales. The fake factor is defined as fhe ratio of fhe number of Thad candidales 
idenlified as medium, fo fhe number salisfying fhe loose, bul nol fhe medium, crileria [61]. The px- 
dependenl fake faclors are measured in dala conlrol samples separately for fhe Thad candidates wilh one or 
Ihree fracks and for IT-i-jels, mulfijef, Z-i-jels and fop quark conlribufions. The fk-i-jels, mulfijef, Z-i-jels and 
lop quark conlrol samples are selecfed by reversing fhe mx cuf (see below), relaxing fhe lepfon isolation 
requiremenf, reversing fhe dileplon velo or by requiring f?-lagged jels, respeclively. The fake faclors 
delermined from Ihese conlrol samples are consislenl wilhin Iheir slalislical uncerlainlies. They are Ihen 
applied fo fhe signal conlrol sample, i.e., evenls passing Ihe selection, excepl lhal Ihe Thad candidate is 


required to satisfy the loose, but not the medium, Thad identification, to estimate the fake Thad background. 
The composition of the sample is determined from a mixture of data-driven methods and MC simulations 
and it is found that the sample is dominated by the IT-i-jets and multijet events. Details of the method 
can be found in Ref. [61]. The method is validated using the same-sign /’Thad events that are otherwise 
selected in the same way as the signal candidates. 

The Z ^ TT background is modeled using selected Z ^ /i/i events in data through embedding [63], where 
the muon tracks and associated energy depositions in the calorimeters are replaced by the corresponding 
simulated signatures of the final-state particles of tau decays. In this approach, the kinematics of the 
produced boson, the hadronic activity of the event (jets and underlying event) as well as pileup are taken 
from data [6]. Other processes passing the Z ^ /r/i selection, primarily from top quark production, are 
subtracted from the embedded data sample using simulation. Their contributions are approximately 2% 
for events with one /^-tagged jet and 25% for events with two or more /^-tagged jets. The Z ^ tt 
background derived is found to be in a good agreement with that obtained from the MC simulation. 

The remaining backgrounds, mostly tt and diboson events with genuine /rhad in their decay final states, 
are estimated using simulation. The small contributions from single SM Higgs boson production and 
from Z{-^ eelfifj.) -i-jets events (in which one of the electrons or muons is misidentified as Thad) are also 
estimated from simulation. The production rates of these processes are normalized to the theoretical cross 
sections discussed in Sec. 2. For the simulation of the tf process, the top quark pj distribution is corrected 
for the observed difference befween data and simulation [64]. The background from misidentified leptons 
is found to be negligible. 

Figures 2(a) and 2(b) compare the observed ditau {nirr) and dijet (mhb) mass distributions with those 
expected from background events after the preselection discussed above. The sample is dominated by 
contributions from top quark production, fake Thad, and Z ^ tt backgrounds. Also shown in the figures are 
the expected signal distributions for a Higgs boson pair production cross section of 20 pb as an illustration. 
The yield of the nonresonant production is significantly higher than that of the resonant production for 
the same cross section, largely due to the harder pT spectrum of the Higgs boson h of the nonresonant 
production. The ditau invariant mass is reconstructed from the electron or muon, Thad, and F™®® using a 
method known as the missing mass calculator (MMC) [65]. The MMC solves an underconstrained system 
of equations with solutions weighted by resolution and the tau-lepton decay topologies. It returns 
the most probable value of the ditau mass, assuming that the observed lepton, Thad and F™'*® stem from a 
TT resonance. The dijet mass is calculated from the two leading f7-tagged jets, or using also the highest-pT 
untagged jet if only one jet is /^-tagged. 

Additional topological requirements are applied to reduce the background. As shown in Fig. 2(c), the 
signal events tend to have small values of the transverse mass nij calculated from the lepton and 
system. Consequently, a requirement of < 60 GeV is applied, which reduces the background signi¬ 
ficantly with only a small loss of the signal efficiency. In addition, the angular separation in the transverse 
plane between the F™'®* and Thad is required to be larger than one radian to reduce the fake Thad back¬ 
ground. 

Background events from tt WWbb —> ivTvbb decay have an identical visible final state to the signal 
if the tau lepton decays hadronically. For signal /j ^ tt —> /Thad events, however, the px of the lepton 
tends to be softer than that of the Thad due to the presence of more neutrinos in the t —> / decays. Thus 
the px of the electron or muon is required to satisfy px(/) < Px(Thad) + 20 GeV. The tt events of the 
tt WW bb —> /v qq' bb final state with a misidentified Thad candidafe remain a large background. To 
reduce ifs contribution, a W boson candidate is reconstructed from the Thad candidate and its closest 
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Figure 2: Kinematic distributions of the hh —> bbrr analysis after the preselection (see text) comparing data with the 
expected background contributions: (a) ditau mass reconstructed using the MMC method, (b) dijet mass mi,/, 
and (c) the transverse mass m^'' of the lepton and the system. The top quark background includes contributions 
from both ff and the single top-quark production. The background category labeled “Others” comprises diboson and 
Z ^ eeliijj. contributions. Contributions from single SM Higgs boson production are included in the background 
estimates, but are too small to be visible on these distributions. As illustrations, expected signal distributions for 
a Higgs boson pair production cross section of 20 pb are overlaid for both nonresonant and resonant Higgs boson 
pair production. A mass of nin - 300 GeV is assumed for the resonant production. The last bin in all distributions 
represents overflows. The gray hatched bands represent the uncertainties on the total background contributions. 
These uncertainties are largely correlated from bin to bin. 
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untagged jet and its mass m-^j is calculated. The W candidate is then paired with a Zj-tagged jet to form 
the top quark candidate with a reconstructed mass nirjb- The pairing is chosen to minimize the mass sum 
m£b + nirb for events with two or more Z^-tagged jets. If only one jet is ^-tagged, one of the ft-jets in the 
mass sum is replaced by the highest-px untagged jet. An elliptical requirement in the form of in the 
Qurj, tfirjb) plane: 

cos 6 - Anit sin 9^ j Am^ sin 6 + Am, cos 0 
28 GeV I 18 GeV 

is applied to reject events with (m^j, nirjb) consistent with the hypothesis (mvy, m,), the masses of the W 
boson and the top quark. Here Amvy = m,-y - mvy, Am, = nirjb - mt and 0 is a rotation angle given by 
tan 6 = nttlmw to take into account the average correlation between nirj and nirjb. 

Finally, the remaining events must have 90 < lUbb <160 GeV, consistent with the expectation for the 
h ^ bb decay. For the nonresonant search, m^r is used as the final discriminant to extract the signal, and 
its distribution is shown in Fig. 3(a). The selection efficiency for the gg ^ hh ^ bbTT signal is 0.57%. 
For the resonant search, the MMC mass is required to be in the range of 100 < m,-,- < 150 GeV. The 
mass resolutions of rubb and m,-,- are comparable for the signal, but the rubb distribution has a longer tail. 
The resonance mass mbbrr reconstructed from the dijet and ditau system is used as the discriminant. To 
improve the mass resolution of the heavy resonances, scale factors of niblmbb and mhlmrr are applied 
respectively to the four-momenta of the bb and tt systems, where nib is set to the value of 125 GeV 
used in the simulation. The resolution of mbbrT is found from simulation studies to vary from 2.4% at 
niH = 260 GeV to 4.8% at 1000 GeV. The improvement in the resolution from the rescaling is largest at 
low mass and varies from approximately a factor of three at 260 GeV to about 30% at 1000 GeV. The 
reconstructed ifibbTr distribution for events passing all the selections is shown in Fig. 3(b). The efficiency 
for the gg^H —>hh—> bbrr signal varies from 0.20% at 260 GeV to 1.5% at 1000 GeV. These efficiencies 
include branching ratios of the tau lepton decays, but not those of heavy or light Higgs bosons. 

To take advantage of different signal-to-background ratios in different kinematic regions, the selected 
events are divided into four categories based on the ditau transverse momentum (less than or greater 
than 100 GeV) and the number of (^-tagged jets (ub = 1 or > 2) for both the nonresonant and resonant 
searches. The numbers of events expected from background processes and observed in the data passing 
the resonant hh^bbrT selection are summarized in Table 2 for each of the four categories. The expected 
number of events from the production of a Higgs boson with mh - 300 GeV and a cross section of 
o-{gg -^H) X BR(H hh) = 1 pb for each category is also shown for comparison. 

Systematic uncertainties from the trigger, luminosity, object identification, background estimate as well 
as Monte Carlo modeling of signal and background processes are taken into account in the background 
estimates and the calculation of signal yields. The impact of these systematic uncertainties varies for 
different background components and event categories. For the most sensitive Ub > 2 categories, the main 
background contributions are from top quark, fake Thad^ and Z —> tt. The jet energy scale and resolution is 
the largest uncertainty for the top-quark contribution, ranging between 10% and 19% for the nonresonant 
and resonant searches. The leading source of systematic uncertainty for the fake Thad background is the 
“fake-factor” determination, due to the uncertainties of the sample composition. Varying the composition 
of IT-i-jets, Z-i-jets, top quark and multijet events in the control samples by ±50% leads to a change 
in the estimated fake T^ad background by 9.5%. The most important source of systematic uncertainty 
for the Z ^ TT background is the ti subtraction from the Z ^ pp sample used for the embedding, due 
to the uncertainty on the tt normalization. Its effect ranges from 8% to 15%. The overall systematic 
uncertainties on the total background contributions to the high (low) category of «& > 2 are 12% 
(9%) for the nonresonant search and 14% (14%) for the resonant search. The largest contributions are 
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Table 2; The numbers of events predicted from background processes and observed in the data passing the final 
selection of the resonant search for the four categories. The top quark background includes contributions from both 
tt and the single top-quark production. The “others” background comprises diboson and Z —> ee/yUyU contributions. 
The numbers of events expected from the production of a nin - 300 GeV Higgs boson with a cross section of 
cr{gg^ H)xBR{Hhh) = 1 pb are also shown as illustrations. The uncertainties shown are the total uncertainties, 
combining statistical and systematic components. 



Hb 

- 1 

Hb 

> 2 

Process 

<100 GeV 

> 100 GeV 

<100 GeV 

> 100 GeV 

SM Higgs 

0.5 ±0.1 

0.8 ±0.1 

0.1 ± 0.1 

0.2 ±0.1 

Top quark 

30.3 ± 3.6 

19.6 ± 2.5 

30.9 ± 3.0 

23.6 ± 2.5 

Z^TT 

38.1 ±4.4 

20.2 ± 3.7 

6.8 ± 1.8 

2.6 ± 1.0 

Fake Thad 

37.0 ± 4.4 

12.1 ± 1.7 

13.7 ± 1.9 

5.4 ± 1.0 

Others 

3.2 ± 3.7 

0.5 ± 1.5 

0.7 ± 1.6 

0.2 ± 0.7 

Total background 

109.1 ± 8.6 

53.1 ±6.0 

52.2 ± 8.2 

32.1 ±5.4 

Data 

92 

46 

35 

35 

Signal rriH = 300 GeV 

0.8 ± 0.2 

0.4 ± 0.2 

1.5 ± 0.3 

0.9 ± 0.2 
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Figure 3: Distributions of the final discriminants used to extract the signal: (a) for the nonresonant search and 
(b) nihhTT for the resonant search. The top quark background includes contributions from both tt and the single 
top-quark production. The background category labeled “Others” comprises diboson and Z —»ee/p/r contributions. 
Contributions from single SM Higgs boson production are included in the background estimates, but are too small 
to be visible on these distributions. As illustrations, the expected signal distributions assume a cross section of 
10 pb for Higgs boson pair production for both the nonresonant and resonant searches. In (b), a resonance mass 
of 300 GeV is assumed. The gray hatched bands represent the uncertainties on the total backgrounds. These 
uncertainties are largely correlated from bin to bin. 
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from jet and tau energy scales and Zj-tagging. The modeling of top quark production is also an important 
systematic uncertainty for the category with two or more Z>-tagged jets and high pl^. 

The uncertainties on the signal acceptances are estimated from experimental as well as theoretical sources. 
The total experimental systematic uncertainties vary between 12% and 24% for the categories with two 
or more Z>-tagged jets, and are dominated by the jet and tau energy scales and Zj-tagging. Theoretical 
uncertainties arise from the choice of parton distribution functions, the renormalization and factorization 
scales as well as the value of strong coupling constant as used to generate the signal events. Uncertainties 
of 3%, 1% and 3% from the three sources, respectively, are assigned to all signal acceptances. 

For the nonresonant search, the observed ditau mass distribution agrees well with that of the estimated 
background events as shown in Fig. 3(a). For the resonant search, a small deficit with a local significance 
of approximafely 2cr is observed in fhe dafa relafive fo fhe background expecfafion af mbbrr ~ 300 GeV 
as is shown in Fig. 3(b). No evidence of Higgs boson pair production is presenf in fhe dafa. The resulting 
upper limifs on Higgs boson pair production from fhese searches are described in Sec. 9. 


7 hh^yyWW* 

This secfion describes fhe search for Higgs boson pair producfion in fhe hh yyWW* decay channel, 
where one Higgs boson decays fo a pair of photons and fhe ofher decays fo a pair of W bosons. The h^yy 
decay is well suifed for lagging fhe Higgs boson. The small Higgs boson widfh fogefher wifh fhe excellenf 
defecfor resolufion for fhe diphofon mass sfrongly suppresses background confribufions. Moreover, fhe 
h WW* decay has fhe largesl branching ratio afler h bb. To reduce mulfijef backgrounds, one of fhe 
W bosons is required to decay fo an elecfron or a muon (eilher direclly or Ihrough a lau leplon) whereas 
fhe ofher is required fo decay hadronically, leading fo fhe yybvqq' final sfafe. 

The dafa used in fhis analysis were recorded wifh diphofon friggers wifh an efficiency close to 100% 
for diphofon evenls passing fhe final offline selection. The diphofon selection follows closely fhaf of fhe 
ATLAS measuremenf of fhe h^yy production rale [56] and lhal of fhe hh —> yybb analysis [21]. Evenls 
are required fo have fwo or more idenlified photons wifh fhe leading and subleading pholon candidales 
having pjlmyy > 0.35 and 0.25, respectively, where niyy is fhe invarianl mass of fhe fwo selecled phofons. 
Only evenls wifh niyy in fhe range of 105 < niyy <160 GeV are considered. 

Addifional requiremenfs are applied fo idenfify fhe h —> WW* —> ivqq' decay signafure. Evenls musl 
have fwo or more jels, and exacfly one lepfon, salisfying fhe idenlificafion crileria described in Sec. 3. 
To reduce mulfijef backgrounds, fhe evenls are required to have wifh significance greafer fhan one. 
Evenls wifh any f?-fagged jel are veloed to reduce confribufions from lop quark production. 

A lofal of 13 evenls pass fhe above selection. The final hh^yyWW* candidales are selecled by requiring 
fhe diphofon mass niyy fo be wifhin a ±2a- window of fhe Higgs boson mass in /j —> 77 where cr is faken fo 
be 1.7 GeV. Due fo fhe small number of evenls, bolh nonresonanf and resonanl searches proceed as counf- 
ing experimenls. The seleclion efficiency for fhe hh yyWW* signal of SM nonresonanf Higgs boson 
pair producfion is eslimaled using simulation to be 2.9%. Eor fhe resonanl producfion, fhe corresponding 
efficiency varies from 1.7% af 260 GeV fo 3.3% al 500 GeV. These efficiencies include fhe branching 
ratios of fhe W boson decays, bul nof fhose of fhe Higgs boson decays. 
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The background contributions considered are single SM Higgs boson production (gluon fusion, vector- 
boson fusion and associated production of Wh, Zh and tth) and continuum backgrounds in the niyy spec¬ 
trum. Events from single Higgs boson production can mimic the hh —> yyWW* signal if, for example, the 
Higgs boson decays to two photons and the rest of the event satisfies the h —> WW* —> ivqq' identification. 
These events would exhibit a diphoton mass peak at m/,. As in the hh bbxT analysis, their contributions 
are estimated from simulation using the SM cross sections [27]. The systematic uncertainty on the total 
yield of these backgrounds is estimated to be 29%, dominated by the modeling of jet production (27%). 
The total number of events expected from single SM Higgs production is therefore 0.25 ± 0.07 with con¬ 
tributions of 0.14, 0.08 and 0.025 events from Wh, tth and Zh processes, respectively. Contributions from 
gluon and vector-boson fusion processes are negligible. 
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Figure 4; The disttibution of the diphoton invariant mass for events passing (a) the relaxed requirements and (b) the 
final selection. The relaxed requirements include all final selections except those on the lepton and F’™^. The red 
curves represent the continuum background contributions and the blue curves include the contributions expected 
from single SM Higgs boson production estimated from simulation. The continuum background contributions in 
the signal niyy mass window are shown as dashed lines. 


The background that is nonresonant in the yy mass spectrum is measured using the continuum background 
in the m^y spectrum. The major source of these backgrounds is Wyy + jets events with a IT —> T’v decay. 
These events are expected to have a diphoton mass distribution with no resonant structure at nih and their 
contribution (N^^) in the signal region, niyy € mh± 2cr, is estimated from the niyy sidebands in the data: 




/sB 

1 - f SB 


Here is the number of events in the data sidebands, defined as the mass region 105 < niyy < 
160 GeV excluding the signal region. The quantity /sb is the fraction of background events in 105 < 
niyy < 160 GeV falling into the signal mass window, and can in principle be determined from a fit of 
the observed niyy distribution to an ansatz function. However, the small number of events after the final 
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selection makes such a fit unsuitable. Instead, /sb is determined in a data control sample, selected as 
the signal sample without the lepton and requirements. Figure 4(a) shows the niyy distribution of 
events in the control sample. For the fit, an exponential function is used to model the sidebands and a 
wider region of nih ± 5 GeV is excluded to minimize potential signal contamination in the sidebands. The 
fit yields a value of /sb = 0.1348 ± 0.0001. Varying the fit range of the sidebands leads to negligible 
changes. Different fit functions, such as a second-order polynomial or an exponentiated second-order 
polynomial, lead to a difference of 1.4% in /sb. To study the sample dependence of /sb, the fit is repeated 
for the control sample without the jet and F"™®® requirements and a difference of only 2% is observed. 
Simulation studies show that the continuum background is dominated by W{€v)yy -i-jets production. The 
■yy^y-i-jets events generated using MadGraph reproduce well the observed niyy distribution. The potential 
difference between yy-i-jets and yyf’y-i-jets samples is studied using simulation. A difference below 1% is 
observed. Taking all these differences as systematic uncertainties, the fraction of background events in the 
signal mass window is /sb = 0.135 ± 0.004. With 9 events observed in the data sidebands, it leads 

to = 1.40 + 0.47 events from the continuum background. Figure 4(a) also shows the contribution 
expected from single SM Higgs boson production. The data prefer a larger cross section than the SM 
prediction for single Higgs boson production, consistent with the measurement reported in Ref. [66]. 

The uncertainties on the signal acceptances are estimated following the same procedure as the hh bbrT 
analysis. The total experimental uncertainty is found to vary between 4% and 7% for different signal 
samples under consideration, dominated by the contribution from the jet energy scale. The theoretical 
uncertainties from PDFs, the renormalization and factorization scales, and the strong coupling constant 
are 3%, 1%, and 3%, respectively, the same as for the hh^bbrT analysis. 

The rtiyy distribution of the selected events in the data is shown in Fig. 4(b). In total, 13 events are 
found with 105 < niyy <160 GeV. Among them, 4 events are in the signal mass window of m/, + 2cr 
compared with 1.65 ± 0.47 events expected from single SM Higgs boson production and continuum 
background processes. The p-value of the background-only hypothesis is 3.8%, corresponding to 1.8 
standard deviations. 

Assuming a cross section of 1 pb (cr{gg hh) or cr{gg H) X BR{H hh)) for Higgs boson pair 
production, the expected number of signal events is 0.64 ± 0.05 for the nonresonant production. For 
the resonant production, the corresponding numbers of events are 0.47 ± 0.05 and 0.72 ± 0.06 for a 
resonance mass of 300 GeV and 500 GeV, respectively. The implications of the search for Higgs boson 
pair production are discussed in Sec. 9. 


8 Combination procedure 


The statistical analysis of the searches is based on the framework described in Refs. [67-70]. Profile- 
likelihood-ratio test statistics are used to measure the compatibility of the background-only hypothesis 
with the observed data and to test the hypothesis of Higgs boson pair production with its cross section 
as the parameter of interest. Additional nuisance parameters are included to take into account systematic 
uncertainties and their correlations. The likelihood is the product of terms of the Poisson probability 
constructed from the final discriminant and of nuisance parameter constraints with either Gaussian, log¬ 
normal, or Poisson distributions. Upper limits on the Higgs boson pair production cross section are 
derived using the CL^ method [71]. For the combinations, systematic uncertainties that affect two or 
more analyses (such as those of luminosity, jet energy scale and resolutions, A-tagging, etc.) are modeled 
with common nuisance parameters. 
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For the hh bbzT analysis, Poisson probability terms are calculated for the four categories separately 
from the mass distributions of the ditau system for the nonresonant search (Fig. 3(a)) and of the bbzT 
system for the resonant search (Fig. 3(b)). The nthtTv distributions of the resonant search are rebinned to 
ensure a sufficient number of events for the background prediction in each bin, in particular a single bin 
is used for mhbrr ^ 400 GeV for each category. For the hh yyWW* analysis, event yields are used 
to calculate Poisson probabilities without exploiting shape information. The hh yybb and hh bbbb 
analyses are published separately in Refs. [21, 22]. However, the results are quoted at slightly different 
values of the Higgs boson mass m/, and, therefore, have been updated using a common mass value of 
mil = 125.4 GeV [24] for the combinations. The decay branching ratios of the Higgs boson h and their 
uncertainties used in the combinations are taken from Ref. [27]. Table 3 is a summary of the number of 
categories and final discriminanfs used for each analysis. 


Table 3: An overview of the number of categories and final discriminant distributions used for both the nonresonant 
and resonant searches. Shown in the last column are the mass ranges of the resonant searches. 


hh 

final sfafe 

Nonresonanf search 

Resonanf search 

Cafegories 

Discriminanf 

Categories 

Discriminanf 

mh [GeV] 

yybb 

1 

TtTyy 

1 

evenf yields 

260-500 

yyWW* 

1 

even! yields 

1 

even! yields 

260-500 

bbTT 

4 

ffljT 

4 

fyibbrr 

260-1000 

bbbb 

1 

even! yields 

1 

^bbbb 

500-1500 


The four individual analyses are sensitive to different kinematic regions of the hh production and decays. 
The combination is performed assuming that the relative contributions of these regions to the total cross 
section are modeled by the MadGraphS [39] program used to simulate the hh production. 


9 Results 

In this section, the limits on the nonresonant and resonant searches are derived. The results of the hh 
bbTT and hh —> yyWW* analyses are first determined and then combined with previously published results 
of the hh —> yybb and hh —> bbbb analyses. The impact of the leading systematic uncertainties is also 
discussed. 

The observed and expected upper limits at 95% CL on the cross section of nonresonant production of 
a Higgs boson pair are shown in Table 4. These limits are to be compared with the SM prediction of 
9.9 ± 1.3 fb [17] for gg^hh production with m/, = 125.4 GeV. Only the gluon fusion production process 
is considered. The observed (expected) cross-section limits are 1.6 (1.3) pb and 11.4 (6.7) pb from the 
hh bbTT and hh yyWW* analyses, respectively. Also shown in the table are the cross-section limits 
relative to the SM expectation. The results are combined with those of the hh —> yybb and hh bbbb 
analyses. The p-value of compatibility of the combination with the SM hypothesis is 4.4%, equivalent to 
1.7 standard deviations. The low p-value is a result of the excess of events observed in the hh yybb 
analysis. The combined observed (expected) upper limit on cr{gg hh) is 0.69 (0.47) pb, corresponding 
to 70 (48) times the cross section predicted by the SM. The hh bbbb analysis has the best expected 
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sensitivity followed by the hh yybb analysis. The observed combined limit is slightly weaker than that 
of the hh^bbbb analysis, largely due to the aforementioned excess. 


Table 4: The expected and observed 95% CL upper limits on the cross sections of nonresonant gg^hh production 
at Vs = 8 TeV from individual analyses and their combinations. SM values are assumed for the h decay branching 
ratios. The cross-section limits normalized to the SM value are also included. 


Analysis 

yybb 

yyWW* 

bbTT 

bbbb 

Combined 

Upper limit on the cross section [pb] 

Expected 

1.0 

6.7 

1.3 

0.62 

0.47 

Observed 

2.2 

11 

1.6 

0.62 

0.69 

Upper limit on the cross section relative to the SM prediction 

Expected 

100 

680 

130 

63 

48 

Observed 

220 

1150 

160 

63 

70 


The impact of systematic uncertainties on the cross-section limits is studied using the signal-strength 
parameter ji, defined as the ratio of the extracted to the assumed signal cross section (times branching ratio 
for the resonant search). The resulting shifts in ju. depend on the actual signal-strength value. 
For illustration, they are evaluated using a cross section of 1 pb for gg —> {H -^hh, comparable to the limits 
set. The effects of the most important uncertainty sources are shown in Table 5. The leading contributions 
are from the background modeling, f7-tagging, the h decay branching ratios, jet and measurements. 
The large impact of the fj-tagging systematic uncertainty reflects the relatively large weight of the hh 
bbbb analysis in the combination. For the hh bbrr analysis alone, the three leading systematic sources 
are the background estimates, jet and F"™®® measurements, and lepton and Thad identifications. For the 
hh yyWW* analysis, they are the background estimates, jet and F'™*® measurements and theoretical 
uncertainties of the decay branching ratios of the Higgs boson h. 

For the resonant production, limits are set on the cross section of gg^H production of the heavy Higgs 
boson times its branching ratio BR(// ^ hh) as a function of the heavy Higgs boson mass mn- The 
observed (expected) limits of the hh bbxT and hh yyWW* analyses are illustrated in Fig. 5 and listed 


Table 5: The impact of the leading systematic uncertainties on the signal-strength parameter yu of a hypothesized 
signal for both the nonresonant and resonant [iriff = 300, 600 GeV) searches. For the signal hypothesis, a Higgs 
boson pair production cross section {cr{gg ^ hh) or cr{gg^H) x BR(i/—>/i/i)) of 1 pb is assumed. 


Nonresonant search 


Resonant search 


Source 

AyU/yU [%] 

niH = 300 GeV 

niH = 600 GeV 

Source 

AyU/yU [%] 

Source 

AyU/yU [%] 

Background model 

11 

Background model 

15 

f^-tagging 

10 

fj-tagging 

7.9 

Jet and 

9.9 

/zBR 

6.3 

/iBR 

5.8 

Lepton and Thad 

6.9 

Jet and 

5.5 

Jet and 

5.5 

ABR 

5.9 

Luminosity 

2.7 

Luminosity 

3.0 

Luminosity 

4.0 

Background model 

2.4 

Total 

16 

Total 

21 

Total 

14 
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in Table 6 (along with results from the hh yybb and hh bbbb analyses). The niH search ranges are 
260-1000 GeV for hh bbrT and 260-500 GeV for hh yyWW*. For the hh bbrr analysis, the 
observed limit around mn ~ 300 GeV is considerably lower than the expectation, reflecting the deficit 
in the observed distribution. At high mass, the limits are correlated since a single bin is used 

for > 400 GeV. The decrease in the limit as mu increases is a direct consequence of increasing 
selection efficiency for the signal. This is also true for the hh —> yylTlT* analysis as the event selection is 
independent of mn- 



300 400 500 600 700 800 900 1000 250 300 350 400 450 500 

[GeV] mu [GeV] 
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Figure 5; The observed and expected upper limit at 95% CL on cr(gg H) x BR(H —> hh) at ^ TeV as a 

function of nin from the resonant (a) hh —> bbrr and (b) hh yyWW* analyses. The search ranges of the resonance 
mass are 260-1000 GeV for hh^ bbrr and 260-500 GeV for hh^ yyWW* . The green and yellow bands represent 
+ lcr and +2cr ranges on the expected limits, respectively. 


Table 6: The expected and observed 95% CL upper limits on crigg—^H) x BR{H ^ hh) in pb at ^fs = TeV from 
individual analyses and their combinations. The SM branching ratios are assumed for the light Higgs boson decay. 


niH Expected limit [pb] Observed limit [pb] 


[GeV] 

yybb 

yyWW* 

bbrr 

bbbb 

Combined 

yybb 

yyWW* 

bbzT 

bbbb 

Combined 

260 

1.70 

11.2 

2.6 

- 

1.1 

2.29 

18.7 

4.2 

- 

2.1 

300 

1.53 

9.3 

3.1 

- 

1.2 

3.54 

15.1 

1.7 


2.0 

350 

1.23 

7.8 

2.2 

_ 

0.89 

1.44 

13.3 

2.8 

- 

1.5 

400 

1.00 

6.9 

0.97 

_ 

0.56 

1.00 

11.5 

1.5 

_ 

0.83 

500 

0.72 

5.9 

0.66 

- 

0.38 

0.71 

10.9 

1.0 

- 

0.61 

500 

— 

— 

0.66 

0.17 

0.16 

— 

— 

1.0 

0.16 

0.18 

600 

- 

- 

0.48 

0.070 

0.067 

- 

- 

0.79 

0.072 

0.079 

700 

_ 

_ 

0.31 

0.041 

0.040 

_ 

_ 

0.61 

0.038 

0.040 

800 

_ 

_ 

0.31 

0.028 

0.028 

_ 

_ 

0.51 

0.046 

0.049 

900 

_ 

_ 

0.30 

0.022 

0.022 

- 

- 

0.48 

0.015 

0.015 

1000 

— 

— 

0.28 

0.018 

0.018 

— 

— 

0.46 

0.011 

0.011 


The hh yybb and hh bbbb analyses are published separately and the mass range covered by the 
two analyses are 260-500 GeV and 500-1500 GeV, respectively. The results of these four analyses, 
summarized in Table 6, are combined for the mass range 260-1000 GeV assuming the SM values of the 
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h decay branching ratios. To reflect the better mass resolutions of the hh bbbb and hh yybb analyses, 
the combination is performed with smaller mass steps than those of the hh bbxT and hh yyWW* 
analyses. The most significant excess in the combined results is at a resonance mass of 300 GeV with 
a local significance of 2.5cr, largely due to the 3.0cr excess observed in the hh —> yybb analysis [21]. 
The upper limit on cr{gg H) x BR{H hh) varies from 2.1 pb at 260 GeV to 0.011 pb at 1000 GeV. 
These limits are shown in Fig. 6 as a function of itih- For the low-mass region of 260-500 GeV, both 
the hh yybb and hh bbrr analyses contribute significantly to the combined sensitivities. Above 
500 GeV, the sensitivity is dominated by the hh^bbbb analysis. Table 5 shows the impact of the leading 
systematic uncertainties for a heavy Higgs boson mass of 300 GeV and 600 GeV. As in the nonresonant 
search, the systematic uncertainties with the largest impact on the sensitivity are from the uncertainties 
on the background modeling, ^-tagging, jet and measurements, and the h decay branching ratios. 
These limits are directly applicable to models such as those of Refs. [72-77] in which the Higgs boson h 
has the same branching ratios as the SM Higgs boson. 



itIh [GeV] 


Figure 6; The observed and expected 95% CL upper limits on cr{gg H) x BR{H hh) at -yT = 8 TeV as 
functions of the heavy Higgs boson mass mn, combining resonant searches in hh yybb, bbbb, bbrr and yyWW* 
final states. The expected limits from individual analyses are also shown. The combination assumes SM values 
for the decay branching ratios of the lighter Higgs boson h. The green and yellow bands represent +lcr and ±2cr 
uncertainty ranges of the expected combined limits. The improvement above nin = 500 GeV is due to the sensitivity 
of the hh —»bbbb analysis. The more finely spaced mass points of the combination reflect the better mass resolutions 
of the hh —> yybb and hh —»bbbb analyses than those of the hh —> bbrr and hh yyWW* analyses. 
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10 Interpretation 


The upper cross-section limits of the resonant search are interpreted in two MSSM scenarios, one referred 
to as the hMSSM [28, 29] and the other as the low-tb-high [30]. In the interpretation, the CP-even light and 
heavy Higgs bosons of the MSSM are assumed to be the Higgs bosons h and H of the search, respectively. 
The natural width of the heavy Higgs boson H where limits are set in these scenarios is sufficiently smaller 
than the experimental resolution, which is at best 1.5%, that its effect can be neglected. 

In the hMSSM scenario, the mass of the light CP-even Higgs boson is fixed to 125 GeV in the whole 
parameter space. This is achieved by implicitly allowing the supersymmetry-breaking scale ms to be very 
large, which is especially true in the low tan p region where ms » 1 TeV, and making assumptions about 
the CP-even Higgs boson mass matrix and its radiative corrections, as well as the Higgs boson coupling 
dependence on the MSSM parameters. Here tan/3 is the ratio of the vacuum expectation values of the 
two doublet Higgs fields. The “low-tb-high” MSSM scenario follows a similar approach, differing in that 
explicit choices are made for the supersymmetry-breaking parameters [30]. The mass of the light Higgs 
boson is not fixed in this scenario, but is approximately 125 GeV in most of the parameter space. The m^ 
value grows gradually from 122 GeV at m^ ~ 220 GeV to 125 GeV as mA approaching infinity. Higgs 
boson production cross sections through the gluon-fusion process are calculated with SusHi 1.4.1 [78-80] 
for both scenarios. Higgs boson decay branching ratios are calculated with HDECAY 6.42 [81] following 
the prescription of Ref. [29] for the hMSSM scenario and with FeynHiggs 2.10.0 [82-84] for the low-tb- 
high scenario. 

The upper limits on (T{gg H)x BR(H hh) can be interpreted as exclusion regions in the (tan/3, 
plane. In both scenarios, the Higgs boson pair production rate (T{gg —> //) x BR(// —> hh) depends on 
tan/3 and the mass of the CP-odd Higgs boson {mA), and so does the mass of the heavy CP-even Higgs 
boson H. The values of and mu are generally different: mu can be as much as 70 GeV above mA in 
the parameter space relevant for this publication with the difference in masses decreasing for increasing 
values of tanyS or mA. Consfant mn lines for a few selected values are shown in Fig. 7. The decay 
branching ratios of the light Higgs boson in these scenarios depend on tan/3 and niA and are different 
from the corresponding SM values used to derive the upper limits shown in Table 6. The upper limits, 
as functions of m//, are recomputed; the hh decay fractions for each final state are fixed to their smallest 
value found in 1 < tan/3 < 2, the range of the expected sensitivity. This approach yields conservative 
limits, but simplifies the computation as the limit calculation does not have to be repeated at each tan p 
value. The results are used to set exclusions in the (tan^S, wta) plane as shown in Fig. 7. The analysis 
is sensitive to the region of low tanyS and mA values in the range ~ 200-350 GeV. For mA < 200 GeV, 
mii is typically below the 2m/, threshold of the H ^ hh decay, whereas above 350 GeV, the H ^ hh 
decay is suppressed because of the dominance of the H ^ tt decay. The observed exclusion region in the 
(tan/3, wta) plane is smaller than the expectation, reflecting the small excess observed in the data. 


11 Summary 


This paper summarizes the search for both nonresonant and resonant Higgs boson pair production in 
proton-proton collisions from approximately 20 fb“^ of data at a center-of-mass energy of 8 TeV recorded 
by the ATFAS detector at the FHC. The search is performed in hh bbxT and yyWW* final states. No 
significant excess is observed in the data beyond the background expectation. Upper limits on the hh 
production cross section are derived. Combining with the hh yybb, bbbb searches, a 95% CF upper 
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Figure 7: The observed and expected 95% CL exclusion regions in the (tan/l, ma) plane of MSSM scenarios from 
the resonant search: (a) the hMSSM scenario and (b) the low-tb-high scenario. The green dotted lines delimit 
the +lcr uncertainty ranges of the expected exclusion regions. The gray dashed lines show the constant values of 
the heavy CP-even Higgs boson mass. The improved sensitivity in the expected exclusion on the contour line of 
niH ~ 260 GeV reflects the improved expected limit on the cross section while the hole or the wedge around the 
niff ~ 325 GeV contour line in the observed exclusion is the result of a small excess at this mass, see Fig. 6. The 
gray shaded region in (b) shows the region where the mass of the light CP-even Higgs boson is inconsistent with 
the measured value of 125.4 GeV. There is no such region in (a) by construction. 


limit of 0.69 pb on the cross section of the nonresonant hh production is observed compared with the 
expected limit of 0.47 pb. This observed upper limit is approximately 70 times the SM gg hh production 
cross section. For the production of a narrow heavy resonance decaying to a pair of light Higgs bosons, 
the observed (expected) upper limit on cr{gg^H) x BR{H^hh) varies from 2.1 (1.1) pb at 260 GeV to 
0.011 (0.018) pb at 1000 GeV. These limits are obtained assuming SM values for the h decay branching 
ratios. Exclusion regions in the parameter space of simplified MSSM scenarios are also derived. 
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Eakultat fiir Physik und Astronomic, Julius-Maximilians-Universitat, Wurzburg, Germany 


42 



Fachbereich C Physik, Bergische Universitat Wuppertal, Wuppertal, Germany 
Department of Physies, Yale University, New Haven CT, United States of Ameriea 
Yerevan Physics Institute, Yerevan, Armenia 

Centre de Calcul de I’lnstitut National de Physique Nucleaire et de Physique des Particules (IN2P3), 
Villeurbanne, France 

“ Also at Department of Physics, King’s College London, London, United Kingdom 
^ Also at Institute of Physics, Azerbaijan Academy of Sciences, Baku, Azerbaijan 
^ Also at Novosibirsk State University, Novosibirsk, Russia 
Also at TRIUMF, Vancouver BC, Canada 

^ Also at Department of Physics, California State University, Fresno CA, United States of America 
^ Also at Department of Physics, University of Fribourg, Fribourg, Switzerland 

^ Also at Departamento de Fisica e Asttonomia, Faculdade de Ciencias, Universidade do Porto, Portugal 
Also at Tomsk State University, Tomsk, Russia 
' Also at CPPM, Aix-Marseille Universite and CNRS/IN2P3, Marseille, France 
■> Also at Universita di Napoli Parthenope, Napoli, Italy 
^ Also at Institute of Particle Physics (IPP), Canada 

^ Also at Particle Physics Department, Rutherford Appleton Laboratory, Didcot, United Kingdom 
Also at Department of Physics, St. Petersburg State Polytechnical University, St. Petersburg, Russia 
” Also at Louisiana Tech University, Ruston LA, United States of America 
° Also at Institucio Catalana de Recerca i Estudis Avancats, ICREA, Barcelona, Spain 
P Also at Department of Physics, The University of Michigan, Ann Arbor MI, United States of America 
^ Also at Graduate School of Science, Osaka University, Osaka, Japan 
'' Also at Department of Physics, National Tsing Hua University, Taiwan 

* Also at Department of Physics, The University of Texas at Austin, Austin TX, United States of 
America 

' Also at Institute of Theoretical Physics, Ilia State University, Tbilisi, Georgia 
“ Also at CERN, Geneva, Switzerland 

" Also at Georgian Technical University (GTU),Tbilisi, Georgia 
Also at Manhattan College, New York NY, United States of America 
^ Also at Hellenic Open University, Patras, Greece 
^ Also at Institute of Physics, Academia Sinica, Taipei, Taiwan 
^ Also at EAE, Universite Paris-Sud and CNRS/IN2P3, Orsay, Prance 

“ Also at Academia Sinica Grid Computing, Institute of Physics, Academia Sinica, Taipei, Taiwan 
Also at School of Physics, Shandong University, Shandong, China 

Also at Moscow Institute of Physics and Technology State University, Dolgoprudny, Russia 
Also at Section de Physique, Universite de Geneve, Geneva, Switzerland 
Also at International School for Advanced Studies (SISSA), Trieste, Italy 

Also at Department of Physics and Astronomy, University of South Carolina, Columbia SC, United 
States of America 

Also at School of Physics and Engineering, Sun Yat-sen University, Guangzhou, China 
Also at Paculty of Physics, M.V.Eomonosov Moscow State University, Moscow, Russia 
Also at National Research Nuclear University MEPhl, Moscow, Russia 
Also at Department of Physics, Stanford University, Stanford CA, United States of America 
Also at Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Budapest, 
Hungary 

Also at University of Malaya, Department of Physics, Kuala Eumpur, Malaysia 

* Deceased 


43 



